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DETERMINATION OF MEAN ATMOSPHERIC DENSITIES 


FROM THE EXPLORER IX SATELLITE 

By Gerald M. Keating^ James A. Mullins^ Claude W. Coffee, 
and David S. McDougal 
Langley Research Center 


SUMMARY 

A method is developed for determining mean atmospheric densities from 
changes in the orbital elements of the Explorer IX (I 96 I Delta l) satellite, a 
12- foot-diameter balloon. The decay of the total energy of the satellite 
orbiting about an oblate spheroid is evaluated and the energy change due to the 
work done by radiation force upon the satellite is subtracted to yield the 
energy decay attributed to atmospheric drag. The direct solar radiation force 
upon Explorer IX is determined from measured vectorial reflectances of the sat- 
ellite material. Equations are given for the effective coefficient of drag of 
Explorer IX, and mean atmospheric densities at heights from 510 km to 7^5 Ion 
are tabulated for the entire lifetime of the satellite (February 
March 1964) from changes in the orbital elements of Explorer IX. 


INTRODUCTION 


L through 


Since the advent of earth satellites, there has been increased interest 
in, and opportunity for, measurement of atmospheric density at extreme alti- 
tudes. Besides the scientific value of using density data to obtain further 
physical insight into the natvire of the atmosphere, the densities must be known 
in order to predict perturbations in the orbits of earth satellites and to 
predict satellite lifetimes. Because of the exceedingly low atmospheric den- 
sities at satellite altitudes, an indirect method is normally employed to infer 
the densities, namely: the density is deduced from the magnitude of atmospheric 

dreig which is computed for the satellite from certain observed orbitel pertur- 
bations. The Explorer IX (196I Delta l) satellite was designed with a low mass- 
area ratio exclusively for the measurement of atmospheric density by the 
indirect method. 

Exactly how the general indirect method is applied depends upon the 
physical characteristics and orbital regime of the satellite. In this paper, 
the specific method used for determination of mean atmospheric densities from 
the Explorer IX satellite is developed and completely described and the den- 
sities thus obtained are tabulated so that Interested parties may analyze them. 
Mean atmospheric densities in the vicinity of the satellite’s perigee are 



deduced from that portion of the energy decay of the satellite which is attrib- 
uted to aerodynamic drag. Expressions and values are determined from the time 
rate of change of energy about an oblate spheroid in terms of mean elements, 
the total radiation force upon the satellite, the energy change due to radia- 
tion force, and the coefficient of drag averaged about one orbit. 

In order to obtain mean atmospheric densities from the energy decay of 
other satellites the satellite mass, average cross section, perturbing forces 
effecting energy decay, thermal accommodation coefficient, satellite wall tem- 
peratures, and many other factors would have to be taken into account. 
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SYMBOLS 


A area 

A-2 = J2^E ~ 66 . 05^6 km 

a semimajor aocis of osculating orbit averaged over 1 revolution 

ao semimajor axis of osculating orbit 

Aaj) change in semimajor axis of orbit per revolution due to aerodynamic 

drag averaged over many revolutions 

Cp) aerodynamic drag coefficient 

Cp effective aerodynamic drag coefficient 

Cs solar constant at 1 astronomical unit 

c speed of light in a vacuum 

D aerodynamic drag 

E energy 

AEp change in energy of satellite per revolution due to aerodynamic drag 

AEp change in energy of satellite per revolution due to aerodynamic drag 

averaged over many revolutions 

AEg energy change of satellite per revolution due to solar radiation 

force 

Bp total energy of satellite 

AEtji total energy change of satellite 
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FS 
H = 

h 

I 

X 

J2 

K 


if 


eccentricity- 

force exerted on satellite by direct solar radiation 




R 




sat 

rO) 

S 


S 


dp/dh 

height above earth 

incident power of light per unit area normal to direction of 
propagation 

orbit inclination referred to equatorial plane of date 

coefficient of second harmonic in earth gravitational potential 

correction factor for rotation of atmosphere 

unit vector in direction of propagation of li^t 

satellite mass 

mean (anomalistic) motion 

coefficients in expression for orbital elements 
generalized orbital elements 
perturbing potential 

distance from center of earth to satellite 
equatorial radius of earth 

distance from earth to sun in astronomical units 
distance from center of earth to point on shadow ellipse 
radius of satellite 

vectorial reflectance of surface for light incident at angle 3 
molecular speed ratio 

molecular speed ratio of diffusely reflected particles 
arc of satellite orbit 
satellite velocity 
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t 

u 

V 

V 

W 

X,Y,Z 

x,y,z 

xi 

a 

OS 

Op 

3 

5S 

5p 

e 

0 

0(3) 

X 

V 

P 

P 

cr 

T 

k 


time 

earth gravitational potential 
potential energy of satellite 
true anomaly 
work 

Cartesian coordinate axes 
Cartesian coordinates 
generalized rectangular coordinate 

change in right ascension between ascending node and perigee point 

right ascension of sun 

right ascension of perigee 

angle of incidence 

declination of sun 

declination of perigee 

eccentric anomaly 

angle measured in orbital plane from semimajor axis of shadow semi- 
ellipse about center of earth in direction of satellite motion 
to re 

angle between normal to area element and direction of vectorial 
reflectance 

angle measured at center of earth from ascending node of satellite 
orbit to earth- sun line 

product of universal gravitational constant and mass of earth 
angle between earth-sun line and its projection on orbital plane 
atmospheric density 

mean atmospheric density 
standard deviation 

angle measured in direction of satellite motion from perigee to pro- 
jection of earth-sun line on orbital plane 


0 

SI 


colatitude angle in spherical coordinate system 
azimuth angle in spherical coordinate system 

right ascension of ascending node referred to equinox of 1950*0 
(i> argument of perigee 

Subscripts: 

in shadow entry 

out shadow exit 

p at perigee 

t at time t 


EXPERIMENT DESCRIPTION 


Because of the exponential decrease of atmospheric density with height, 
most of the change of energy due to atmospheric drag occurs in the vicinity of 
perigee for moderate orbital eccentricities. The higher the eccentricity of 
the satellite orbit, the greater is the portion of energy decay in the immedi- 
ate vicinity of perigee, and, therefore, the more acciirately the position of the 
density measurement may be pinpointed. The lower the mass-area ratio of the 
satellite, the greater is the variation of the orbit due to energy decay and, 
therefore, the more accurate the density may be measured. 

In order to obtain accurate values for both the magnitude and position of 
a density determination, it was, therefore, advantageous to place a satellite 
with a low mass-area ratio into an eccentric orbit. The Explorer IX satellite 
(fig. l) was designed specifically for this purpose. The I2. 00-foot-diameter 
satellite weired Ik. 62 pounds, had a lower mass-area ratio than all previous 
satellites except Echo I (i960 Iota), and yet was sturdy enough to remain spher- 
ical at heights above I50 km (ref. l). Because it was spherical, the satellite 
presented the same known frontal area to the atmosphere regardless of its 
orientation. Also it coiild be tracked optically and by radar because of its 
shiny metallic surface. Furthermore, the lifetime of the experiment is not 
limited to electronic Instrument lifetime but to satellite lifetime, because 
the entire satellite served as a sensor. 

Explorer IX was launched on February I6, I96I, from the NASA Wallops 
Station by means of the solid-propellant Scout vehicle. The initial orbit had 
an inclination of 38 * 8 ®^ a hel^t of perigee of 652 km, and a height of apogee 
of 2586 km. It was predicted that such an orbit would be close to a resonant 
condition for radiation force which covild substantially reduce the satellite 
lifetime. Therefore, the satellite was launched in the early morning, allowing 
the perigee altitude to increase for 1 year and to decrease for the next 
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Figure 1.- Explorer IX satellite. L-61-5T62 


2 years. If an evening launch had been attempted^ radiation force would have 
driven the satellite into the earth in less than a year. Besides increasing 
the satellite lifetime^ a morning launch made it possible to obtain densities 
for the same altitude at two different times. 


PROCEDURE 


The procedure used to obtain mean atmospheric densities may be outlined 
as follows: 

(1) Accurate observations of the satellite were continually obtained from 
various sources. 

(2) Time-dependent orbital elements were calculated from these observa- 
tions by means of a differential improvement program. 

(3) The energy decay rate was deduced from the time-dependent orbital 
elements. 
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(4) By eliminating the computed energy changes due to forces other than 
atmospheric drag^ the energy decay rate due to drag was deduced. 

( 5 ) A mean atmospheric density was deduced from the energy decay rate due 
to drag and certain aerodynamic properties of the satellite when certain 
properties of the atmosphere were assimied. 

These steps are discussed in greater detail in the following sections. 


Obsei^ations 

The shiny metallic surface of the balloon permits both radar and optical 
tracking. The Baker- Nunn camera network and moonwatch teams were the primary 
sources of position measurements used to determine orbital elements. Publica- 
tions are available which describe the Baker-Nunn cameras of the Smithsonian 
Institution Astrophysical Observatory (for example^ ref. 2) and the methods of 
reduction of the photographs to position measurements (for example, ref. 
therefore, only a brief discussion is given herein. 

For bright satellites such as Explorer IX a fixed camera with a rotating 
barrel shutter is used. This shutter interrupts the trace of the satellite 
five times during an exposure (fig. 2). At the instant at which the central 
break occurs, the time is recorded on the film. Satellite direction relative 
to the star background at the time of the photograph is determined by measuring 
the center of the central break with respect to background star images. These 
observations are accurate to a few minutes of arc when field- reduced, and they 
are later reduced with accuracies of a few seconds of arc. (See ref. 5*) An 
average of about 1 observation per 
satellite revolution was obtained 
for the lifetime of the satellite. 

Another important source of obser- 
vations is the numerous moon- 
watch teams throughout the world. 

Moonwatch observations, though 
accurate to only about 0.5^^ are 
an important supplement to the 
Baker-Nunn data. 


Determination of Orbits 

The Smithsonian Astrophysi- 
cal Observatory Differential Orbit 
Improvement Program (refs, k 
and 5 ) '^as used to determine the 
orbital elements of the satellite 
from the observations. The method 
of determination of the orbital 
elements is essentially as 
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follows. A set of expressions for the orbital elements is assumed. These 
elements vary with time and are composed of polynomials with the option of 
additional transcendental and exponential terms^ such as 




where 


qj orbital elements 

constants 

t time from epoch 

It is these constants Pj ^ which are to be improved. 

The expressions for the elements are evaluated at times corresponding to 
the times at which observations were made and short-term gravitational pertur- 
bations (ref. 6) are added. The observations are transformed to rectangular 
sidereal coordinates and the computed satellite positions are compared with 
the observed positions. The residual from a given measurement can be equated 
with a linear combination of corrections for the elements . The partial deriva- 
tives of the rectangular coordinates with respect to the elements are coeffi- 
cients in the following resultant relation: 




rectangular coordinate 


(Axj_)^ residual in coordinate x^ at time t 


The increments Aqj can in turn be written as linear combinations of the error 
in the coefficients follows: 



If enough measurements are made, the corrections determined. In 

practice, many more observations are made than the minimum number necessary, 
and the corrections ^j,k determined by a method of least squares. 

The standard deviation cr of the residuals is then computed and all meas- 
urements with residuals greater than 3cr are discarded. The entire process is 
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repeated until all the residuals are within the acceptable maximum. At 
this time the program prints out a set of parameters 

In practice^ least -mean- square fits devoid of short-periodic perturba- 
tions are found for eccentricity e, orbital inclination i, argument of 
perigee cu, right ascension of ascending node 9 ., and mean anomaly M. The 
semimajor axis a is determined from the following relation (ref. 6): 



where n^ the mean (anomalistic) motion, is dM/dt. The elements derived by 
the above procedure are referred to as mean elements. 


Satellite Energy in Terms of Mean Elements 

A relationship between the energy of the satellite and the mean elements 
can be obtained in the following manner. The total energy of the satellite 

is the sum of the kinetic energy ^ ms'^ and the potential energy V. The 

kinetic energy for the oscvilating ellipse is given by ^(— - ) and the 

2 \r aoy 

potential energy due to the gravitational potential of the earth U is -Urn, 
where 


and 


Then 



1 p 1 ? 

i sin'^i + — sin i cos 


2(o> + v) 



i sin^i + i sin% 


cos 2(co + v) 



or 




( 5 ) 


( 6 ) 


Bp 



( 7 ) 


Equation ( 7 ) gives the energy for the osculating ellipse in terms of the 
osculating semimajor axis and the perturbing potential. If the satellite is 
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perturbed only by the oblateness potential, the energy for the osculating 
ellipse is constant and therefore equal to the energy for the mean ellipse. 
According to reference 6, the semimajor axis of the mean ellipse is related to 
the semimajor axis of the osculating ellipse by 



( 8 ) 


Substitution of equation ( 8 ) and equation (5) into equation (7) yields 


pn 


- 1 sln%){l - ' 

2a 

L 3a^ ' 

2 /'■ ' 


(9) 


Energy changes in the orbit can be deduced from changes in the orbital elements. 
Since the energy change is observed to be small. 


AErj. = ^ Aa + 
2el 


^(1 - 


sin i cos i Ai 



^ sin%^(l - e^) 



( 10 ) 

In the orbital regime of Explorer IX, the changes of e and i are small 
enough that the energy change may be approximated by the equation 

AEt « ^ Aa (11) 

2a^ 

with errors of the order of J2 or about 0.1 percent. Therefore, a change of 
energy may be measured in terms of a change in the mean semimajor axis. 


Causes of Energy Change 

In order to evaluate the energy change due solely to atmospheric drag, 
the energy changes due to other perturbing forces must be determined and sub- 
tracted from the total energy change. 

Gravitational effects .- If the satellite is perturbed only by the earth 
gravitational field, the satellite total energy with respect to the earth 
remains constant. The lunar and solar gravitational perturbations may be 
ass\imed not to affect the semimajor axis (refs. 7 and 8). Therefore, from 
equation (ll) the change in energy due to these effects is negligible. 

Charge drag .- Because of the high thermal velocity of electrons relative 
to satellite velocity and the lower ion velocity, the satellite acquires a 
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negative charge which is decreased to some extent by photoemission. Conse- 
quently, ions which would not have been affected by a neutral satellite strike 
or are deflected by the satellite. The resulting drag, additional to that 
which would have been experienced by a neutral satellite, is referred to as 
electrostatic drag which has been estimated for the Explorer IX satellite 
(ref. 9) 'to be about 2 percent of the drag upon a neutral satellite at I 5 OO km 
for average solar activity with no photoemission assumed. Photoemission reduces 
the charge and therefore the relative effect of electrostatic drag. It should 
be noted that the percentage ionization decreases with decreasing solar activity 
and decreasing altitude. Since the perigee point of Explorer IX was continu- 
ously in the sunlight throioghout its lifetime, and the solar activity and alti- 
tude were less than those assumed, it is believed that this effect would be 
less than 1 percent. 

The satellite becomes polarized when cutting the earth’s magnetic lines 
of force. With an increased electron flux on the positively charged side of 
the satellite and an increased ion flux on the negatively charged side, a cur- 
rent flows across the satellite which interacts with the earth's magnetic field 
to produce an induction drag. At I 5 OO km (ref. S), the induction drag has been 
determined to be about 1 percent of the neutral drag. For the range of heights 
of the Explorer IX perigee, this effect is quite negligible. Therefore, the 
change of energy due to electrostatic and electromagnetic drag forces will lead 
to errors of less than 1 percent in atmospheric density. 

Radiation forces .- Finally, the effects of radiation forces should be con- 
sidered. The change of energy produced by direct solar radiation has been 
found to be of the same order as drag on Explorer IX and is much larger than 
the changes due to earth- reflected radiation, earth-emitted radiation, and 
satellite-emitted radiation. Besides being relatively small in magnitude, 
earth- reflected radiation and earth-emitted radiation produce forces upon the 
Explorer IX satellite which are essentially radial. Because of the thermal 
design of the satellite, the temperature of surface elements facing the sun is 
approximately the same as the temperature of surface elements on the opposite 
side of the satellite. (See ref. 10.) Therefore, the net force produced by 
satellite-emitted radiation is quite small. 

The energy change in the satellite due to direct solar radiation depends 
upon the motion of the satellite with respect to the earth shadow and upon the 

force Fg exerted by solar radiation. The force acting on flat samples of 
Explorer IX material due to the incidence and reflection of radiation may be 
determined from reference 11 and Fg can be determined by integration over the 
satellite surface. 

From reference 11, the force dE^g on a flat surface due to the incidence 
and reflection of solar radiation is 

.l^[S.?(|3)]dA (12) 
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where 


I incident power per unit area normal to direction of 

propagation 

3 angle between Incident ray and normal to surface 

c speed of light in a vacuum, taken as 2.997 95 k 10^^ cm/sec 

unit vector in direction of incident light 
"r(p) vectorial reflectance of surface at particular angle of incidence 

Vectorial reflectance: A particle of electromagnetic radiation has asso- 

ciated with it momentum which is transferred to an absorbing or reflecting 
mediiim. If parallel light of energy E is incident on a surface at rate 
dE/dt the magnitude of the force exerted on the surface due to incidence as 

opposed to reflection is ^ ^ "the direction of the incident light. The 

resultant force exerted on the surface by reflected light is dependent upon the 
amount of light departing in each direction. The vector siunmation of the rate 
of momentum exchange by the reflected light divided by the magnitude of the 
rate of momentum exchange by the incident light is vectorial reflectance. In 
other words, vectorial reflectance is the ratio of the vectorial summation of 
reflected power to the magnitude of the parallel incident power. 

In reference 11, the vectorial reflectance of the Explorer IX material was 
evaluated by using simulated sunlight, 
the angle at which the light was inci- 
dent. The results given in refer- 
ence 11, together with equation (12), 
can be used to obtain the total force 
on the spherical satellite. 

Total direct solar radiation 
force: Let parallel light of inten- 

sity I illiuninate a hemisphere 
shown in figure 5^ be an incre- 

mental surface area element of the 
hemisphere, and the radius of the 
hemisphere which is parallel to the 
incident light define the Z-axis of a 
Cartesian coordinate system origi- 
nating at the center of the hemi- 
sphere. The orientation of the X- 
and Y-axes is arbitrary. Points on 
the hemisphere can be described by 
the spherical coordinates 0 and 
\|f, where 0 is measured from the 
Z-axis and \[r is measured from the 
X-axis to the plane containing 0. 


The vectorial reflectance depended upon 


Z 



Figure Geometry of momentum transfer by 
solar radiation. 
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The force on the hemisphere can be found by resolving the force on each 
surface element into components normal to and parallel with the Z-axis and 
integrating these forces over the hemisphere. From simple geometrical con- 
siderations it can be shown that ^ equals <ji, the resultant force due to_^ 
symmetry is directed negatively along the Z-axis, and the z- component of r(p) 
is 


|cos [0 + 0(p^ if (13) 

Here 0(3), the angle between the normal to dA. and 'r(3)> is assumed to be 
in the plane of 0. The area element for a sphere of radius ^ga,t becomes 

dA = sin 0 d0 d\jf (l^^) 

Proper substitution of relation (13) s^nd equation (lU) into equation (12) 
and integrating over the hemisphere yields 



The Explorer IX is unifonnly covered with 2- inch-diameter dots of white paint 
which constitute I8 percent of its svirface area. In equation (16) 

cos [0 + 0(0)] be replaced by the weighted vectorial reflectance 
conponents 


O.82{jl?(0)|cos[0 + e(0)]}^^inum P + ®^*^^)white paint 

The vectorial reflectances given in reference 11 with extrapolations to 0° 
and 90° incidence were substituted into the above expression. A curve of the 
integrand in equation (16) was then grajdiically integrated to yield 

(17) 
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If the incident light is sunlight, 

I = ^ (18) 

"■I 

where 

Cg solar constant at 1 astronomical unit, taken as 1.395 X 10^ erg/cm^-sec 

rg distance from earth to sun in astronomical vinits 

Energy change per revolution; The change in^ energy of the satellite jier 
revolution due to direct solar radiation force Fg may be given as follows: 



where ds is an increment of arc along the satellite orbit. Let the Z-axis 
of a geocentric coordinate system be the projection of the earth- sun line on 
the orbital plane (fig. 4). Let the angle between the earth- sun line and the 
Z-axis be v and the angle between the perigee point on the orbit and the 
Z-axis measured in the direction of satellite motion be t. Because of the 
slow variation of the direction of the sun relative to the osculating orbit 
during 1 revolution t and v are assumed constant. Due to the minute varia- 
tion in the intensity of the solar flux |^^sj assumed constant during 

1 revolution when the satellite is in sunlight. If the satellite enters the 
earth shadow at and exits shadow at ZQut^ follows that 

AEg = |fs|cos 

( 20 ) 

where and are 

determined from ecjua- 
tions ( 35 ) and ( 36 ). 

Although the satel- 
lite actually passes 
through the penxunbra 
before going into complete 
shadow, for the sake of 
calculation, a geometrical 
shadow is assumed about a 
spherical earth. The axis 
of the cylindrical shadow 
is colinear with the earth- 
sun line. The intersec- 
tion of the orbital plane 



Figure 4.- Shadow and orbit ellipses. 
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with the cylindrical shadow of the earth is a semiellipse. The intersection 
of this shadow semiellipse with the orbit ellipse determines the orbital posi- 
tions of shadow entry and exit. 


The shadow semiellipse, centered at the center of the earth, has semiminor 

Ttp 

axis r-p and semimajor axis of — - — . The radius of an ellipse with center 
^ sin V 

at the origin can be written in polar coordinates as 


r2 = 


a^b^ 

a^sin^e + b^cos^G 


( 21 ) 


where a and 
radius of the 
from 


b are semimajor and semiminor axes, respectively. Then the 
shadow semiellipse (denoted by subscript e) may be determined 



_E 

2 - 2 P 

sin 0 + sin V cos 9 


( 22 ) 


where -90° 9 ^ 90°. The angle 9 is measured from the semimajor axis of 

the semiellipse to re in the direction of satellite motion. 


For the purposes of obtaining intersection points, the orbit of the satel- 
lite can be adequately expressed as an ellipse with the center of the earth at 
one focus as follows: 


a(l - e^) 

r = — ^ — 

1 + e cos V 


(25) 


where 


' V = 9 + T + 180° (24) 

I The angles v and t may be determined from considerations of the geo- 

I metrical relationship between the orbital plane, the sun, and the perigee point. 

I Shown in figure 5 is the celestial sphere upon which are shown the intersection 

I of the orbital plane, the intersection of lines passing from the center of the 
earth throxigh the sun S, the ascending node AN, the perigee P, the normal 
to the orbital plane J., and the projection of the earth-sun line on the 
orbital plane Z. The vernal equinox VE and north celestial pole NP are 
also illustrated. 

Referring to spherical triangle 

cos(90° - v) = cos i cos ^90° - + 

■ Therefore, 

V = sin”^|cos i sin 8g 

i 

I 15 

I 

I 


S-NP-J., from the law of cosines 

sin i sin ^90° - (°S ~ ^ + 90°) (25) 

- sin i cos 5g sin^og - (26) 


I 


where 


-90° V ^ 90° (27) 

Equation (26) may be substituted into equation (22) to obtain as a func- 

tion of 0 and the known quantities. 

Referring again to spherical triangle S-KP-1, from the law of cosines 

cos (90° ~ ^s) “ cos(90° - v)cos i + sin(90° - v)sin i cos (90° - t - go) (28) 

therefore. 


sin(r + go) 


sin 63 - sin v cos i 
sin i cos v 


(29) 


In order to determine the quadrant of t + go, an additional relationship may be 
developed for cos(t + go). From spherical triangle S-J.-AN, 

cos A = cos(90° - v)cos 90° + sin(90° - v)sin 90° cos(t + go) (50) 

therefore, 

cos(t + go) = ^ (51) 

cos V 

From spherical triangle S-NP-AN, 

cos X = cos ^90° - 6g^cos 90° + sln^90° - Sg^sin 90° cos (52) 

Simplifying equation (52) and substituting into equation (3I) gives 

cos(t + go) 

cos 5s cos(o3 - 
cos V 

(33) 

From equation (27), cos v 
is always positive. Since 
cos 63 is always positive, 

the sign of cos^og - ft) 

is always the sign of 
cos(t + go). Knowing the 
signs of cos(t + go) and 
sin(r + go), the quadrant 
of (t + go) and, therefore, 
the value of (t + go), are 
uniquely determined. 
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From equation { 29 ), 


, / sin 6q - sin v cos i \ 

T = sin-1 S__ , _ ^ 

\ sin 1 cos V / 

Equations (5^), (24)^ and (2J) may be used to obtain r as a function of 0 
and the known quantities. 

Now r and rg are known functions of 9. The values of 9 for which 

r = rg were found by an iterative solution of a quartic equation. Comparisons 
of r and rg were made systematically at different values of 9 until the 

angles at which r equals rg were found within 0.01°. If the satellite 
enters shadow at exits shadow at ^Q^^^Ogut; then 


z 


in 


^out 


Tin ®in 

(55) 

’out ®out 

( 56 ) 


Energy change between epochs: Equations have now been developed for 

determining the energy change due to direct solar radiation during 1 revolu- 
tion. In order to evaluate the energy change due to solar radiation between 
successive sets of orbital elements, the energy change per revolution is eval- 
uated once daily and the total energy change is determined by numerical inte- 
gration. Values of a, e, i, cn, W, og, and 5g used to determine 

AEs are assumed to vary linearly between epochs. This assumption introduces 
a very small error, since o), fl, and otg (the most sensitive values) do, for 

all practical purposes, vary linearly over a time interval of the order of a 
week. If the satellite enters continuous sunlight, causing the energy change 
due to solar radiation to approach zero, or if the energy change per revolution 
changes sign, the energy change per revolution during that day is evaluated 
every 2 hours. 


Calculation of Densities 

If aerodynamic drag D is acting on a satellite traveling along arc s, 
the work W done per revolution upon the satellite by drag may be given as 

w = ^ D ds (37) 

The work done is equivalent to the change of energy per revolution pro- 
duced by drag AEj). The drag is ^s ACp where p is the atmospheric density, 
s is the satellite velocity, A is the cross-sectional area, and Cp is the 
drag coefficient. The energy change AEj) is first determined in terms of true 
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anomaly and then put in terms of eccentric anomaly. If Cp and A are 
assumed constant 


AEd = - 


“2“ 



(1 + e cos 
(1 - e cos 


d€ 


(58) 


vhere p is the mass of the earth times the universal gravitational constant 
taken as 5.986 155 2?! x 10^0 dyne s cm e is eccentricity, and e is 

eccentric anomaly. 

An equation has been derived in reference 12 for the change of semimajor 
axis due to drag in a spherically symmetric atmosphere from an equation similar 
to equation (58) with the exception that a correction factor K has been 
applied to yield the satellite velocity relative to an atmosphere rotating with 
the earth. 


Assiaming that 


-r+r-r 


p = Pp exp 


H 


(59) 


where H is a constant and the subscript p denotes values at perigee, and 
expanding the Integrand in equation (58) as a power series in eccentricity, the 
expression may be replaced by Bessel functions which when in turn are replaced 
by their asymptotic expansions yield the following relation (see ref. 12): 


5 - -1 H^D|.|2e 


p 2 KCd(AA a 


1 _ 2e + ^ - T^-fl - lOe + -?2_) + . 

Ibae 


m 


where 


and 


e< 0.2 

a 


(i^i) 


K = n - ^ cos il 


(42) 


Here, cn is the angular velocity of the earth and Sp can be obtained from 


•2 _ /2 1 ' 

" “ "b ■ ", 


(>•3) 


The density at the perigee point of the satellite may be determined once 
the values of Aap, m/A, H, a, e, Cp, and K are found. 
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The change of semimajor axis per revolution due to drag averaged between 
successive sets of orbital elements is The average change in energy per 

revolution due to drag AEp is given as the difference between the total 

energy change and the energy change due to direct solar radiation between suc- 
cessive se ts o f orbital elements divided by the number of revolutions. The 
value of Aap is then determined from the following equation: 



The greater the nvtmber of revolutions, the less is the uncertainty in 
and, therefore, in The only disadvantage is that due to the preces- 

sion of perigee and of the ascending node, the position at which the density is 
measured becomes more uncertain with an increasing number of revolutions. 
Therefore, a time interval should be chosen such that the density is of suit- 
able accuracy without undue sacrifice of the accuracy of time and position of 
the measurement. 

It has been shown (ref. 13 ) that for a sphere in free molecular flow, the 
drag coefficient for diffuse reflection is 


Cp = 


2 exp(-Sg) /^ 
V 





(^ 5 ) 


The molecular speed ratio S is the ratio of the satellite speed with 
respect to the atmosphere to the most probable molecular speed, and Sr is 
the molecular speed ratio of the diffusely reflected particles. The satellite 
speed depends upon the orbit and the position of the satellite in that orbit. 
The rotation of the atmosphere must be considered in determining satellite 
speed with respect to the atmosphere. The atmospheric mean molecular speed at 
any position depends upon the temperature and molecular weight, which is deter- 
mined by the heating and heat transfer in the atmosphere. The mean molecular 
speed of the diffusely reflected particles depends in addition upon the accom- 
modation coefficient and temperature of the satellite. This means that Cp 
is, under some Initial assumptions for atmospheric conditions and properties 
of the satellite, a function of satellite position, and velocity, which are in 
turn functions of the orbital elements. That is, Cp should not be assvimed 

constant and removed from the integral of equation ( 38 ). 

Because of the difficulty of integrating in equation (38) with Cp a 

function of a1 1 its im plicit variables, an investigation was made with the 
object of determining Cp which could be written in terms of fewer variables. 
Let Cp be defined by 
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AEj) = - 



C„p dC = - 

(l - e cos e) ' 


pACjj 

2 


L 


€=2n 

€=0 


(l + e cos , 

i i—j- pde 

(l - e cos 


(i^6) 


The effects of orbital parameters, thermal accommodation coefficient, satel- 
lite skin temperature, rotating atmosphere, and solar heating were investigated 
by choosing several sets of initial conditions and calculating Cp for a 

variety of orbits. It was found that the most significant variables were 
height of perigee and thermal accommodation coefficient. Assuming skin tem- 
peratures obtained when Explorer IX is in continuous sunlight (ref. lO) and 
complete thermal accommodation, the following approximate equations were deter- 
mined giving Cj) as a function of height of perigee hpi 

^ = 2.050 + 2.484 X 10"^ hp (200 km < hp 6OO km) (4?) 

^ = 2.154 + 1.125 X 10-^ hp (600 km ^ hp < 800 km) (48) 


The value of H, which is the negative slope of a log-density — height 
curve, was obtained from the densities given in the I962 U.S. Standard 
Atmosphere ( ref . 14 ) . 

H = (-5.99^ X 10"5)h2 + 0. 1659b + 7.1687 (200 km ^ h ^ 8OO km) (49) 


Because H is variable, a weighted value of H was used. The 

3Hp 

calculated at the height hp + where Hp is the value of 

(See ref. I5.) The height of perigee in kilometers is given as 


quantity H is 
H at hp . 


hp = a(l - e) - 6378.388(1 - 


sin^io’ 


297 


(50) 


The position of the perigee point, which is essentially the density point, 
may be given in terms of the radius, right ascension, and declination of peri- 
gee as follows: 


rp = a(l - e) 

(51) 

cxp = n + a 

(52) 


where a is the change in right ascension between the ascending node and the 
perigee point. 
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It can be shown by spherical trigonometry that 


a = tan“^(cos i tan co) ( 55 ) 

If 0 ^ i < 90°^ a, and o) are in the same quadrant. If 90 ° ^ i < l 80 °^ a 
is in the quadrant of -tn. 

Similarly, it is found that 


where 


6p = sin"^(sin i sin cu) 

-90° ^ 6p ^ 90° 


( 54 ) 

( 55 ) 


ATMOSPHERIC DENSITIES 


The equations given in the section "Procedure" have been incorporated in 
a computing program for determination of atmospheric density from the energy 
decay of a satellite. Given two sets of orbital elements and the corresponding 
right ascension and declination of the sun, and the earth- sun distance in 
astronomical units, the program determines the following quantities; 

(1) Total energy change due to solar radiation force between epochs 

(2) Total energy change between epochs 

(5) Average energy change per revolution due to drag and the corresponding 
change of semimajor axis per revolution due to drag 

( 4 ) Aerodynamic drag coefficient averaged over a revolution 

(5) The quantity H at h-n + — 

^ 4 

( 6 ) Average values of the orbital elements, the right ascension, declina- 
tion, radius, and height above the geoid of the density point, and parameters 
associated with analysis of density 

(7) Average density and log2_Q of average density 

The general organization of this program is shown in figure 6 . 

Presented in table I are mean orbital elements computed by the Smithsonian 
Astrophysical Observatory using the Differential Orbit Improvement Program 
discussed previously. The elements are referenced to the equinox of 1950*0 and 
the equator of date. The elements tabialated for I96I and 1962 are given in 
references I6, I7, and I8 and are based on precisely reduced Baker-Nunn photo- 
graphs. At this time, precise elements are published for only I96I and I962. 
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Figure 6.- Schematic diagram of computational procedure. 

The elements for dates between January 1^ 19^5^ September 30^ 19^3^ 
given in references 19 ; 20^ and 21 and are based upon less accurate field- 
reduced data. The remainder of the tabulated orbital elements are unpublished 
predictions which were originally used by the Smithsonian Astrophysical 
Observatory to predict satellite position so that subseq^uent observations could 
be made. Consequently^ all these observations were obtained prior to the epoch 
times . 

Given in table l(a) are the published orbital elements of Explorer IX, In 
the first column^ the epoch indicates the reference time for a set of orbital 
elements. All times are zero Universal Time to the nearest millisecond. The 
second column gives the mean (anomalistic) motion in revolutions per day. The 
last four columns give the eccentricity; orbital inclination in degrees^ argu- 
ment of perigee in degrees^ and right ascension of the ascending node in 
degrees. With the exception of the epoch^ the last place in each column is 
generally unreliable. Table l(b) is similar to table l(a) except that the 
second column is semimajor axis in kilometers. 

With the equations given in this paper^ mean atmospheric densities and 
their associated time and position were computed between the successive epochs 
tabulated in table I. Table II gives the results of these computations. The 
first four columns are simply midpoints in time^ radius, right ascension, and 
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declination of the perigee point between successive epochs. The fifth column 
gives log2_o of mean atmospheric density. In the sixth column is the time 
interval between successive orbital elements used to infer atmospheric density. 

Referring to equation (^O), errors in the determination of average density 
may be attributed chiefly to inaccuracies in Aap, Cp, K, and the 

atmospheric model. Any of these factors could introduce an error on the order 
of 5 percent. The error in Aa-Q is principally a function of the accuracy of 

orbit determination^ and the accuracies of determination of energy changes by 
causes other than drag. The errors in orbit determination are dependent upon 
the number of observations, the accuracy of observations, the distribution of 
observations in space and time, and inadequacies in theory. The errors due to 
energy changes other than drag are principally dependent on errors in determi- 
nation of perturbing forces acting on the satellite. Errors in Cj) are prin- 
cipally dependent upon the inaccuracy in the estimate of the thermal accommoda- 
tion coefficient. Errors in are dependent upon errors in the mean 

molecular mass and kinetic temperature. Errors in K depend principally upon 
errors in the rate of atmospheric rotation. Errors due to the atmospheric 
model could result from the choice of a radially symmetric atmosphere instead 
of one with a diurnal bulge. 


CONCLUDING REMARKS 


A specific procedure has been developed and successfully used for deter- 
mination of mean atmospheric densities from changes in the orbital elements of 
the Explorer IX satellite. Errors in determination of mean density may be 
attributed chiefly to inaccuracies in the determination of energy decay due to 
drag, aerodynamic drag coefficient, variation of density with altitude, the 
rate of rotation of the atmosphere, and the assumed atmospheric model. During 
the development of the method, equations were obtained for the total energy of 
a satellite about an oblate spheroid; the solar radiation force and resultant 
energy changes upon Explorer IX; and the aerodynamic drag coefficient of 
Explorer IX as a function of height of perigee. From orbital elements derived 
from observations of Explorer IX, mean atmospheric densities and their asso- 
ciated position were deduced for the entire lifetime of the satellite. 


Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., March l 6 , 1965* 
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Epoch 

n, rev/day 

e 

i, deg 

0), deg 

Ti, deg 

196 1 
FEB . 

18.0 

12.159994 

. 122108 

38.862 

. 106.541 

165.730 

FEB. 

2A.C 

12.16C292 

. 121341 

38. 862 

134.998 

143.890 

MARCH 

2 . C 

12 . 16C617 

. 12C48C 

38. 866 

163.558 

122.063 

MARC )- 

8 .C 

12.161C10 

. 1 19687 

38.869 

192.165 

100.242 

MARCH 

14. C 

12 .161297 

. 118986 

38.867 

220.821 

78.428 

MARCH 

20.0 

12.161522 

; .118344 

38.866 

249.468 

56.615 

MARCH 

26. C 

12 . 161944 

. 1 1 7738 

38.863 

27B. 144 

34.807 

APRIL 

l.C 

12.162977 

1 .117234 

38.862 

306.795* 

13.000 

APRIL 

7.C 

12.163538 

. 1 16881 

38. 860 

335.509 

351.196 

APRIL 

13. C 

12.163632 

. 1 16659 

36.856 

4.159 

329.395 

APRIL 

19.0 

12.164225 

. 1 1 6456 

38.853 

32.734 

307.596 

APRIL 

23. C 

12.164546 

. 116222 

j 38. 849 

61.232 

285.799 

MAY 

l.O 

12.165CC6 

. 115817 

1 38.846 

89.729 

263.995 

MAY 

7.C 

12.1654C7 

. 1 1 5283 

38.846 

11B.218 

242.198 

MAY 

13. C 

12.165598 

. 114733 

38.849 

146.780 

220.407 

MAY 

19. C 

12.165919 

.114211 

38. 85C 

175.391 

198.621 

MAY 

25.0 

12.166781 

. 1 13759 

36. 650 

1 204.072 

176.830 

MAY 

31. C 

12.167636 

. 113316 

38.846 

232. 757 

155.036 

JUNE 

6.0 ' 

12.1664C3 

. 1 12661 

38.847 

261.417 

133.240 

JUNE 

12. C 

12.166695 

. 112425 

38.845 

290.068 

111.444 

JUNE 

Iti .C 

12. 16672 6 

. 112C47 

33.843 

318.732 

89.657 

JUNE 

24. C 

12.169191 

.111744 

38.841 

347.356 

67.877 

JUNE 

28.0 

12.1692C2 

.111653 

38.84C 

6.465 

53.354 

JULY 

4.C 

12. 16911C 

. 111515 

38.839 

35.052 

31.578 

JULY 

10. C 

12.169129 

.111363 

38.835 

63.571 

9.804 

JULY 

16.0 

12 . 16958 1 

. 11 1CC3 

38.834 

92.077 

348.027 

JULY 

22.0 

12.17C159 

. 11C578 

38.835 

120.579 

326.255 

JULY 

28. C 

12.17C7C6 

.11C1C8 

38.837 

149.160 

304.488 

AUG. 

3.C 

12. 17C96C 

. 1C9654 

33.838 

177.801 

282.721 

AUG. 

9.0 

12. 17133C 

• IC9268 

38.838 

206.488 

260.953 

AUG. 

15. C 

12. 172233 

. IC8896 

38.837 

235.216 

239.179 

AUG. 

21.0 

12.173167 

.108466 

38.836 

263. 958 

217.405 

AUG. 

27.0 

12.17339 1 

. 108115 

38.834 

292.686 

195.631 

SEPT. 

2.0 

12.173522 

. 1C7673 

38.831 

321.419 

173.863 

SEPT. 

8.0 

12.172663 

.1C77S0 

38.829 

350.141 

152.097 

SEPT. 

14. C 

12.174249 

. 1C7788 

38.827 

18.841 

130.333 

SEPT. 

20.0 

12. 174912 

. 1C7846 

38.826 

47.453 

108.567 

SEPT. 

26. C 

12.175326 

. 107798 

38.823 

76.001 

86.797 

CCT. 

2.0 

12.17564C 

. 1C75C8 

38.822 

104.509 

65.027 

OCT. 

8.C 

12.176224 

. 1C7C64 

38.823 

133.048 

43.257 

OCT. 

14.0 

12.176684 

. 106558 

38.825 

161.669 

21.492 

OCT. 

2C.C 

12 .1772 72 

. 1C6C81 

38.828 

190.377 

359.725 

OCT. 

26. C 

12.177646 

. 1C5699 

38.828 

219.148 

337.962 

NOV. 

1.0 

12.17839C ; 

. 105418 

38.830 

247.968 

316.201 

NOV. 

7. C 

12.176744 1 

. 105241 

38.328 

276.807 

294.437 

NOV . 

13. C 

12.179341 : 

. 1C5172 

38. 827 

305.640 

272.671 

NOV. 

19.0 

12.179656 1 

. 1C527C 

38.825 

334.458 

250.902 

NOV. 

25. C 

12.16C497 

. IC5536 

38. 822 

3.239 

229.132 

CEC. 

1.0 

12.1813C5 1 

. 1C59C9 

38.821 

31.929 

207.354 

DEC. 

7.C 

12.162361 i 

. 1C6221 

38.816 

60.527 

185.568 

CEC. 

13. C 

12.183131 

. 106259 

38.815 

89.019 

163.770 

CEC. 

19. C 

12.183616 

. 1C5978 

38. 816 

117.552 

141.969 

CEC . 

25.0 

12. 10391 6 

• 105498 

38. 820 

146.130 

120.179 

CEC. 

3 1.0 

12.184163 

. IC4956 

38.821 

174.803 

98.393 












Epoch 

1962 


JAN. 

6.C 

JAN. 

12.0 

JAN. 

18.0 

JAN. 

24. C 

JAN. 

30. C 

FEB. 

5.C 

FEB. 

11. C 

FEB. 

17. C 

FEB. 

23. C 

MARCH 

l.C 

MARCH 

7.C 

MARCH 

13. C 

MARCH 

19. C 

MARCH 

25.0 

MARCH 

31. C 

APRIL 

6.0 

APRIL 

12.0 

APRIL 

18.0 

APRIL 

24.0 

APRIL 

30. C 

MAY 

6.C 

MAY 

12.0 

MAY 

18.0 

MAY 

24.0 

MAY 

3C.C 

JUNE 

5.C 

JUNE 

11.0 

JUNE 

17. C 

JUNE 

23.0 

JUNE 

29.0 

JULY 

5.0 

JULY 

11.0 

JULY 

17. C 

JULY 

23.0 

JULY 

29. C 

AUG. 

4.C 

AUG. 

10.0 

AUG. 

16. C 

AUG . 

22. C 

AUG. 

28.0 

SEPT . 

3.C 

SEPT. 

9.0 

SEPT. 

15.0 

SEPT. 

21. C 

SEPT. 

27. C 

GCT. 

3.C 

OCT. 

9.0 

OCT. 

15. C 

OCT. 

21.0 

OCT. 

27. C 

NOV. 

2.C 

NCV. 

8.C 

NOV. 

14.0 

NOV. 

20. C 

NOV. 

26.0 

EEC. 

2.0 

DEC. 

8.0 

DEC. 

14. C 

DEC . 

20.0 

DEC . 

26. C 


TABLE I.- ORBITAL ELEMEOTS OF EXPLORER IX - Continued 
(a) With mean motion - Continued 


n, rev/day 


12 .164244 
12.184313 
12 .1845CS 
12. 18462S 
12.1853C1 
12.185822 
l2.iee5Cl 
12.187274 
12. ieec7s 
12.188897 
12.189551 
12.189927 
12. 19C462 
12.191575 
12.192739 
12.193438 
12.194129 
12.195C59 
12. 196353 
12.197255 
12.1980CC 
12,198781 
12.199658 
12.20C623 
12.2C16C4 
12.202724 
12.203866 
12.2C4788 
12,205369 
12.205948 
12.2064CC 
12.2C6992 
12.2C7645 
12.2C847C 
12.2C9356 
12.21C176 
12.211087 
12.212C66 
12.213376 
12,214674 
12.215943 
15< 
12.218449 
12.219882 
12.221432 
12.22334C 
12.22518C 
12.227C89 
12.228650 
4i: 
12.232321 
12.2342C8 
12. 2361CC 
12.238181 
12,240252 
12.242195 
12.244377 
12.246682 
12.246947 
12.25C787 


. 1C4532 
.104206 
. 1C4C55 
. 1C3999 
. 1C4C52 
. 104287 
. 1C4625 
.104996 
.1C5225 
. IC5212 
. 1C4976 
. 1C465C 
.104391 
.104239 
. 104248 
. 1C4341 
. 1C4456 
. 1C4613 
. 1C4839 
. 1C5128 
. 1C5478 
. 1C5696 
. 1C5804 
.1C5780 
. 105727 
.105792 
.105921 
.106179 
.106431 
. 1C6639 
.1C686C 
. 107148 
. 1C7524 
.107871 
. 108127 
.108265 
.108278 
. 108332 
. 1C8476 
.1C67C8 
.1C9C61 
.109460 
. 1C9914 
. 11C415 
. 111025 
.111613 
. 112123 
. 112460 
. 1126C4 
. 112607 
. 112628 
. 112764 
.113090 
. 113510 
. 114060 
. 114766 
. 115494 
.116243 
. 116891 
.117359 


38. 623 
38. S25 

38.828 

38.826 

38.824 

33.822 

38. 820 

38. 819 

38.817 

33.817 

38.817 

38.818 

38.820 

38.822 

38.825 

38.826 

38.828 

38.826 

36.825 

38.821 

38.820 

38.820 

38.820 

38.820 

38.822 

38.825 

38.828 

38.828 

38.829 

38.830 
38.832 
38.829 
38.628 

38.827 

38.828 
38.330 
38.832 

38.835 

38.836 
38.838 
38.841 
38.343 
38.843 
38.841 
38.838 

38.837 
38.338 
38. 342 
38.848 
38.853 

38.857 
38. 357 

38.858 
38. 357 
38.360 
38.857 
38. 853 
38.351 
38. 852 
38. 855 







TABLE I.- ORBITAL ELEMENTS OF EXPLORER IX - Continued 
(a) With mean motion - Concluded 


Epoch 

n, rev/day 

e 

1, deg 

0), deg 

rt, deg 

1S63 

JAN. l.C 

12.2516CC 

.117510 

38. 855 

133.940 

201.421 

JAN. 7.C 

12.2525CC 

.117640 

38.866 

163.040 

179.246 

JAN. 13. C 

12.2536EC 

. 11761C 

38. 871 

192. 110 

157.042 

JAN. 19.0 

12.25551C 

.117670 

38.873 

221.292 

134.854 

JAN. 25.0 

12.25726C 

.117970 

38.876 

250.540 

112.654 

JAN. 31. C 

12.2591CC 

.118200 

38. 875 

279.800 

90.440 

FEB. 6.C 

12.2613<^C 

.118950 

38. 882 

3C8 .990 

68.215 

FEB. 12.0 

12.26377C 

.119630 

38.883 

338- 183 

A5.986 

FEB. 18. C 

12.26626C 

.120210 

38. 879 

7.350 

23.730 

FEB. 24.0 

12.26867C 

. 121110 

38.872 

36.480 

1.456 

MARCH 2.0 

12.271 IOC 

.12174C 

38.873 

65.541 

339.167 

MARCH 8.0 

12.27357C 

. 122310 

38.875 

94.570 

316.854 

MARCH 14.0 

12.27651C 

. 122610 

38.875 

123.581 

294.544 

MARCH 2C.C 

12.27926C 

.122830 

38.886 

152.678 

272. 190 

MARCH 26.0 

12.28211C 

.122950 

38.885 

181.876 

249.839 

APRIL l.C 

12.285C1C 

. 123180 

38. 889 

211.163 

227.476 

APRIL 7.C 

12.28816C 

.123480 

38.893 

240. 5C2 

205.087 

APRIL 13. C 

12 .29124C 

.123970 

38.894 

269.834 

182.691 

APRIL IS.C 

12 .29A6CA 

. 124380 

38. 900 

299.147 

160.268 

APRIL 26. C 

12.2S743C 

. 124940 

38. 901 

328.453 

137.837 

MAY l.C 

12.2SS96C 

.12560C 

38.905 

357.790 

115.382 

MAY 7.C 

12 .3C4C6C 

. 126110 

38.903 

27.110 

92.915 

MAY 13.0 

12.3C932C 

. 126640 

38.903 

56.424 

70.432 

MAY 19.0 

12.31413C 

. 127020 

38.902 

85.717 

47.911 

MAY 25. C 

12.317960 

. 1274C0 

38.900 

115.020 

25.379 

MAY 31.0 

12.32219C 

. 127660 

38.901 

144.341 

2.811 

JUNE 6.C 

12.326630 

. 12781C 

38.912 

173.760 

340.214 

JUNE 12.0 

12.33236C 

. 127850 

38.909 

203.268 

317.605 

JUNE 18.0 

12 .33745C 

. 128170 

38.912 

232.818 

294.966 

JUNE 24.0 

12.34144C 

. 128410 

38.915 

262.365 

272.309 

JUNE 30.0 

12.34492C 

.128980 

38.918 

291 .990 

249.623 

JULY 6..C 

12.347987 

. 129370 

38.922 

321.510 

226.915 

JULY 12. C 

12.3517S2 

.129780 

38. 923 

351. 115 

204.195 

JULY 18.0 

12.355935 

. 130320 

38.925 

20.752 

181.467 

JULY 24. C 

12.361639 

. 130830 

38. 925 

50.360 

158.714 

JULY 30.0 

12.36787C 

. 131210 

38.919 

79.983 

135.915 

AUG. 5.0 

12.375603 

. 131550 

38.924 

109.580 

113. 103 

AUG. ll.O 

12.382150 

. 131200 

38.925 

139.160 

90.240 

AUG. 17. C 

12. 3878CC 

. 131500 

38.919 

168.830 

67.370 

AUG. 23. C 

12. 39614C 

. 131490 

38.930 

198.794 

44.436 

AUG. 29.0 

12.4C917C 

. 131530 

38.928 

228.760 

21.450 

SEPT. 4.C 

12.41779C 

. 131 700 

38.938 

258.750 

358.460 

SEPT. 10. C 

12.42556C 

. 131700 

38.945 

288.700 

335.430 

SEPT. 16. C 

12 .436 760 

. 132150 

38.942 

318.850 

312.320 

SEPT. 22.0 

12.45490C 

.132200 

38.923 

349.200 

289.210 

SEPT. 26. C 

12.47353C 

.1316C0 ; 

1 

38.930 

9.100 

273.700 




TABLE I.- ORBITAL ELEMENTS OF EXPLORER IX - Concluded 


(b) With semimajor axis 


Epoch 

km 

e 

1, deg 

(Xi, deg 

deg 

1963 

SEPT. 

26.0 

7852.5669 

.131600 

38.930 

9.100 

273.700 

OCT. 

1.0 

7845.4729 

.131490 

38.951 

34.158 

254.291 

CCT. 

8.0 

7839.6069 

.132010 

38.956 

69.209 

227.111 

OCT. 

15.0 

7828.7819 

.133560 

38.938 

IOA.667 

199.698 

OCT. 

19.0 

7818.8660 

.130886 

38.939 

125.258 

183.917 

OCT. 

26.0 

7804.1389 

.129706 

38.939 

161.112 

156.319 

NOV. 

2.0 

7789.9409 

.128173 

38.957 

197.033 

128.528 

NOV. 

9.0 

7773.7079 

.126771 

38.992 

233.383 

100.585 

NCV. 

16.0 

7756.7869 

.125452 

38.964 

270.099 

72.443 

NCV. 

24.0 

7738.1760 

. 124398 

38.960 

312.325 

40.046 

DEC. 

1.0 

7720.4029 

.124019 

38.960 

349.687 

11.462 

DEC. 

7.0 

7703.2549 

.122791 

38.958 

21.544 

346.795 

DEC. 

12.0 

7689.8359 

.122757 

39.003 

A8.588 

326.098 

DEC. 

17.0 

7676.5609 

. 120772 

38.948 

75.161 

305.296 

DEC. 

23.0 

7661.6589 

.119657 

38.947 

107.474 

280.175 

DEC. 

30.0 

7645.4559 

.117921 

38.953 

145.502 

250.684 

196A 






JAN. 

3.0 

7634.3359 | 

.116690 1 

38.954 

167.455 

233.744 

JAN. 

11.0 

7610.6199 

.114210 

38.964 

211.668 

199.597 

JAN. 

17.0 

7591.1629 

.112310 

38.963 

245.339 

173.761 

JAN. 

25.0 

7565.9399 

.110030 

38.967 

290.692 

139.007 

FEB. 

4.0 

7532.8420 

.107810 

38.963 

347.947 

94.982 

FEB. 

10.0 

7505.0839 

. 104990 

38.951 

22.422 

68.279 

FEB. 

23.0 

7433.5900 

.098400 

38.947 

98.948 

9.066 

MARCH 

1.0 

7376.7449 

.091926 

38.950 

141.046 

336.387 

MARCH 

8.0 

7312.5769 

.083690 

38.949 

184.404 

302.730 

MARCH 

14.0 

7256.7180 

.077040 

38.959 

222.812 

273.125 

MARCH 

22.0 

7182.5720 

.068716 

38.925 

276.095 

232.484 

MARCH 

29.0 

7090.5769 

.057570 

38.940 

324.202 

195.609 


.29 



TABLE II.- MEAIJ DENSITIES AND THEIR ASSOCIATED TIME AND POSITION OF MEASUREMENT AND 
TIME INTERVAL OVER WHICH DENSITIES ARB AVERAGED 

Epoch Op, deg 6p, deg rp, km logpQ (p^ g/cn^)| 


FEB, 

21. C 

282. 213 

32.625 

7014.742 

-16.39 

6.0 

M4RCH 

-l.O 

288. 146 

18.698 

7021.124 

-16.41 

6.C 

MARCH 

5.C 

289.488 

1. 344 

7027.591 

-16.43 

6.0 

MARCH 

11.0 

290.545 

-16.254 

7033.426 

-16.51 

6.0 

MARCH 

17. C 

295. 7C9 

-3C. 991 

7038.690 

-16.52 

6.0 

MARCH 

23. C 

3C7. 775 

-38.595 

7043.549 

-16.42 

6.0 

APRIL 

-2.0 

321.878 

-35.439 

7047.700 

-16.25 

6.0 

APRIL 

A.O 

330.005 

-23. 1 79 

7050.814 

-16.38 

6.0 

APRIL 

10.0 

332.342 

-6. 361 

7052.945 

-16.51 

6.0 

APRIL 

16. C 

333.056 

11. 449 

7054. 509 

-16.49 

6.0 

/tPRIL 

2 . C 

336.548 

27.300 

7056. 115 

-16.52 

6.0 

MAY 

-2.C 

346. 5C3 

37.388 

7058. 515 

-16.46 

6.0 

MAY 

A.C 

0. ece 

37.494 

7062.097 

-16. 5C 

6.0 

MAY 

10.0 

10.933 

27.547 

7066.311 

-16.66 

6.0 

MAY 

16. C 

14.567 

1 1. 734 

7070.492 

-16.69 

6.0 

MAY 

22.0 

15.328 

-6. 085 

7074.151 

-16.56 

6.0 

MAY 

28.0 

17.627 

-22.937 

7077. 391 

-16.61 

6.0 

JUNE 

3.0 

25.642 

-35.291 

7080.661 

-16.66 

6.0 

JUNE 

9.C 

39. 692 

-38.615 

7084.012 

-16.69 

6.0 

JUNE 

15.0 

5 1. 864 

-31. 168 

7087. 198 

-16.65 

6.0 

JUNE 

21.0 

57. 153 

-16. 52C 

7089. 82C 

-16.46 

6.0 

JUNE 

26. C 

58.2C4 

-1. 941 

7091 . 301 

-16.69 

4.0 

JULY 

l.C 

58.914 

12. 843 

7092.231 

-16.82 

6.0 

JULY 

7.C 

62. 867 

28.392 

7093. 402 

-16.77 

6.0 

JULY 

13.0 

73. 427 

37.804 

7095.354 

-16.55 

6.0 

JULY 

IS.C 

87. 745 

36.998 

7098. 287 

-16.51 

6.0 

JULY 

25. C 

97. 322 

26.388 

7101.641 

-16.57 

6.0 

AUG. 

-0. 

100.591 

10.273 

7105. 173 

-16.8? 

6.C 

AUG. 

6.C 

101.347 

-7. 579 

7108.404 

- 16.84 

6.0 

AUG. 

12. C 

1C4.C24 

-24.215 

7111. 132 

-16. 6C 

6.0 

AUG. 

18.0 

112.748 

-35.994 

7114.017 

-16.56 

6.0 

AUG. 

2A.C 

127. 146 

-38.352 

7116.908 

-16.75 

6.0 

SEPT. 

-1.0 

138.845 

-30. C31 

7119.213 

-16.72 

6.0 

SEPT. 

5.0 

143. 663 

-14.909 

7120.453 

-16.75 

6.0 

SEPT. 

11. C 

144.716 

2.314 

7120.650 

-16.60 

6.0 

SEPT. 

17. C 

146.415 

2C.C48 

7120. 187 

-16.60 

6.0 

SEPT. 

23. C 

153.062 

33.514 

. 7119.936 

-16.75 : 

6.0 

OCT. 

-l.O 

166.239 

38.822 

7121.104 

-16.65 

6.0 

OCT. 

5.C 

179.326 

33. 332 

7123. 776 

-16.66 

6.0 

OCT. 

1 l.O 

165. 854 

19. 766 

7127.363 

-16.54 

6.0 

OCT. 

17.0 

187.509 

2.494 

7131.163 

-16.71 

6.C 

OCT. 

23. C 

168 .6C4 

-15.223 

7134.441 

-16.76 

6.0 

NOV. 

-2.C 

193. 6C9 

-3C. 289 

7136.868 

-16.59 

6.0 

NOV. 

A.C 

2C5. 577 

-38.423 

7138.431 

-16.83 

6.0 

NOV. 

IC.C 

220.043 

-35. 766 

7139.276 

- 16 .68 

6.0 

NOV. 

16. C 

228.653 

-23. 743 

7138.943 

-16.79 

6.0 

NOV. 

22.0 

231.281 

-6. 968 

7137.265 

-16.80 

6.0 

CEC. 

-2.C 

232. 112 

IC.917 

7134.433 

-16.75 

6.0 

CEC. 

A.O 

235. 562 

26. 914 

7131. 332 

-16.54 

6.3 

CEC. 

IC.C 

245.412 

37. 215 

7129. 580 

-16.66 

6.0 

CEC. 

16. C 

259. 728 

37.593 

7130. 307 

-16.37 

6.0 

DEC. 

22 .C 

270. C45 

27.041 

7133. 190 

-16.85 

6.0 

DEC . 

28 .0 

273. 835 

12.101 

7137. 160 

-16.76 

6.0 




TABLE II,- MEAN DENSITIES AND THEIR ASSOCIATED TIME AND POSITION OF MEASUREMENT AND 
TIME INTERVAL OVER WHICH DENSITIES ARE AVERAGED - Continued 



Op, deg 

5p/ <ieg 

rp, km 

loglO (p> g/cn>5) 



1S62 

JAN. 3.C 
JAN. 9.C 
JAN. 15.C 
JAN. ^1.0 
JAN. 27, C 


F £0 • 

. 2 • C 

FES. 

8.0 

FEB. 

14.0 

FE6. 

20.0 

MARCH 

-2.0 

MARCH 

4.C 

MARCH 

10.0 

MARCH 

16. C 

MARCH 

22.0 

MARCH 

28.0 

APRIL 

3.0 

APRIL 

9.C 

APRIL 

15. C 

APRIL 

21. C 

APRIL 

27.0 

MAY 

3.0 

MAY 

9.C 

MAY 

15.0 

MAY 

21. C 

MAY 

27. C 

JUNE 

2.C 

JUNE 

8.0 

JUNE 

14. C 

JUNE 

20.0 

JUNE 

26.0 

JULY 

2.0 

JULY 

8.0 

JULY 

14. C 

JULY 

20.0 

JULY 

26.0 

AUG. 

1.0 

AUG. 

7.0 

AUG. 

13.0 

AUG. 

19.0 

AUG. 

25. C 

SEPT. 

-0. 

SEPT. 

6.0 

SEPT. 

12.0 

SEPT. 

18.0 

SEPT. 

24.0 

OCT. 

-0. 

OCT. 

6.0 

OCT. 

12.0 

OCT. 

18.0 

CCT. 

24.C 

NOV. 

-l.C 

NOV. 

5.0 

NOV. 

11.0 

NOV. 

17.0 


DEC. -l.C 
DEC. 5.C 
DEC. il.O 
DEC. 17.0 
CEC. 23.0 


-5.742 
-22.7CC 
-35.206 
-3b. 596 
-21.062 
-16. 265 
1.425 
18.835 
32.745 
38. 795 
34.000 
2C. 782 
3.593 
-14.228 
-29. 578 
-38.237 
-36. 173 
-24.503 
-7. 842 
10.098 
26.307 
36.985 
37.765 
28.233 
12.538 
-5.337 
-22.365 
-35.010 
-38.649 
-31.349 
-16. 644 
1.066 
18.571 
32.621 
38. 800 
34. 059 
20. 818 
3.582 
-14,282 
-29.656 
-38.277 
-36. 118 
-24.337 
-7.571 
10.437 
26.636 
37. 159 
37.645 
27. 841 
11.949 
-6.059 
-23. 108 
-35.514 
-38.524 
-30.511 
-15.345 
2.565 
20.004 
33.607 
38.850 


-16.99 
-17.16 
-17.09 
-16.97 
-I 6.S5 
-i6.ee 
-16.83 
-16.80 
-16.71 
-16.65 
-16.71 
-16.90 
-16.86 
-16.59 
-16.60 
-16.83 
-16.74 
-16.60 
-16.46 
-16.67 
-16.84 
-16.86 
-16,80 
-16.75 
-16.79 
-16.77 
-16.78 
-16.78 
-16.78 
-16.81 
-16.91 
-16.85 
-16.92 
-16,09 
-16.90 
-16.99 
-16.91 
-16.88 
-16.72 
-16.75 
-16.73 
-16.65 
-16.58 
-16.56 
-16.57 
-16.48 
-16. 52 
-16.45 
-16.42 
-16.35 
-16.39 
-16,44 
-16.49 
-16.44 

-16.45 

-16.51 

-16.44 

-16.42 

-16.43 

-16.53 
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TABLE II.- MEAN DENSITIES AND THEIR ASSOCIATED TIME ATH) POSITION OF MEASUREMENT AND 
TIME INTERVAL OVER WHICH DENSITIES ARE AVERAGED - Continued 


Op, deg 


loglO (p^ g/cm3) 


33.C79 
19. UA 
1.52 3 
-16.377 
-31. 317 
-38. 711 
-3A. 866 
-21.680 
-4.537 
13.547 
29. 196 
38. 183 
36.382 
24.770 
7.949 
-10.281 
-26. 765 
-37. 369 
-37.442 
-26.955 
-10.509 
7. 781 
24. 728 
36.444 
38. 150 
28.902 
12.978 
-5.334 
-22.770 
-35. 500 
-38. 555 
-30. 226 
-14.623 
3.724 
21.429 
34. 763 
38.761 
31.239 
15.991 
-2.392 
-20.445 
-34.3C6 
-38. 844 
-31. 500 
-15.979 
-0.538 


>16.68 
-16. 57 
-16.48 
-16.41 
-16.54 
-16.56 
-16.40 
-16.36 
-16.33 
-16.37 
-16.38 
-16.37 
-16.23 
-16.24 
-16.21 
-16.22 
-16.19 
-16.22 
-16.15 
-16.21 
-16.15 
-15.93 
-15.85 
-15.92 
-16.06 
-16.01 
-15.95 
-15.85 
-15.90 
-16.04 
-16.11 
-16.15 
-15.93 
-15.86 
-15.75 
-15.78 
-15.66 
-15.75 
-15.83 
-15.51 
-15.48 
-15.61 
-15.66 
-15.47 
-15.24 
-15.04 


TABLE II.- MEAN DENSITIES AND THEIR ASSOCIATED TIME AND POSITION OF MEASUREMENT AND 
TIME INTERVAL OVER WHICH DENSITIES ARE AVERAGED - Concluded 


Epoch 

Op, deg 

6p, deg 

rp, km 

logic (p^ g/cm5) 

Time interval, 
days 

IS63 

OCT. -1.5 

281.136 

13.395 

6816.521 

-15.18 

5.0 

CCT. ^.5 

285.242 

29. 555 

6809.285 

-15.43 

7.0 

OCT. 11.5 

2SS.465 

38.881 

67<53.931 

-15.14 

7.0 

OCT. 17. C 

312. 7C7 

34.735 

6789.334 

-14.92 

4.0 

OCT. 22.5 

31S.S1C 

22.125 

6793.695 

-15.00 

7.0 

NOV. -1.5 

321. 7C3 

C.585 

6791.693 

-15.01 

7.0 

NOV. 5.5 

323. 3C5 

-21.26C 

6789.860 

-14.95 

7.0 

NOV. 12.5 

333.517 

-36.679 

6785.960 

-14.94 

7.0 

NOV. 2C.C 

352. 76S 

-35. 889 

6779.627 

-14.96 

8.0 

CEC. -2.5 

2.437 

-17.748 

6769. 2A6 

-14.91 

7.0 

CEC. A.O 

3.455 

3.527 

6760. 151 

-14.85 

6.0 

CEC. ^.5 

5.C6C 

21.186 

6751.610 

-14.88 

5.0 

CEC. 1^.5 

ii.iec 

33.691 

6747.658 

-14.88 

5.0 

CEC. 20.0 

24.423 

38.935 

6747.172 

-14.91 

6.0 

CEC. 2E.5 

38.S66 

30. 360 

6744.399 

-14.94 

7.0 

1S6^ 

JAN. 1.0 

43.506 

14. 533 

6743.694 

-14.86 

4.0 

JAN. 7.C 

44.127 

-5.993 

6742.463 

-14.83 

8.0 

JAN. 14. C 

47.987 

-28.C97 

6740.014 

-14.80 

6.0 

JAN. 21. C 

63.825 

-38.937 

6736.044 

-14.81 

8.0 

fEB. -l.C 

83.233 

-24.202 

6727.111 

-14.79 

10.0 

FEB. 7.C 

85.666 

3.257 

6718.945 

-14.64 

6.0 

FEB. 16.5 

92.843 

33.238 

67C9.743 

-14.57 

13.0 

MARCH -2.5 

119. 3C6 

32.986 

6700.469 

-14.40 

7.0 

MARCH 4.5 

125.957 

10. 761 

6699.741 

-14.35 

7.0 

MARCH 11.0 

126.694 

-14. 531 

6699.217 

-14.36 

6.0 

MARCH 18. C 

137. C69 

-36.053 

6693.492 

-14.38 

8.0 

MARCH 25.5 

16C.787 

-32.917 

6685.950 

-14.25 

7.0 


NASA -Langley, 1965 L-45^5 
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